de Haas-van Alphen effect investigation of the 
electronic structure of Al substituted MgB 2 



in 
o 
o 



S3 
O 
O 

Or 
C3 



^3 

C 

o 

(N 
> 

o 
in 
o 

-I— > 



c 

o 
o 



X 



A. Carrington, a J.D. Fletcher, J.R. Cooper, 6 O.J. Taylor, L. Balicas, c 
N.D. Zhigadlo, d S.M. Kazakov, d J. Karpinski, d J.P.H. Charmant 6 and J. Kortus / 
a H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, BS8 1TL, United Kingdom. 
b Department of Physics and IRC in Superconductivity, University of Cambridge, 
Madingley Road, Cambridge CBS ORE, United Kingdom. 
c National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, U.S.A. 
d Laboratorium fur Festkorperphysik, ETH Zurich, CH-8093 Zurich, Switzerland. 
e School of Chemistry, University of Bristol, Bristol BS8 ITS, United Kingdom, and 
f Institut de Physique et Chemie de Materiaux de Strasbourg, Strasbourg, France. 

(Dated: February 2, 2008) 

We report a de Haas- van Alphen (dHvA) study of the electronic structure of Al doped crystals of 
MgB2. We have measured crystals with ~ 7.5% Al which have a T c of 33.6 K, (~ 14% lower than 
pure MgEb). dHvA frequencies for the a tube orbits in the doped samples are lower than in pure 
MgB2, implying a 16 ± 2% reduction in the number of holes in this sheet of Fermi surface. The 
mass of the quasiparticles on the larger a orbit is lighter than the pure case indicating a reduction 
in electron-phonon coupling constant A. These observations are compared with band structure 
calculations, and found to be in excellent agreement. 



The binary compound MgB2 is an unusual supercon- 
ductor exhibiting a very high transition temperature 
(T c — 39 K) and two distinct superconducting gaps. This 
behavior is believed to originate from a phonon mediated 
mechanism with different coupling strengths to the elec- 
tronic a and ir bands m£ The two-gap structure is only 
observable experimentally because of very weak scatter- 
ing between the a and ir bands£ 

Atomic substitutions change the properties of MgB2 
both by increasing the scattering rates and by changing 
the electron density. This provides a way to engineer 
its superconducting properties and also allows us to test 
theoretical understanding of this material. The two el- 
ements which substitute most readily in MgB2 are Al 
and C which replace Mg and B respectively giving the 
general formula Mgi_ a: Al x (Bi_yCj / )2. Both dopants add 
electrons to the material, and cause T c to decrease in a 
similar way^ (i.e., dT c /dx ~ \dT c jdy). 

In principle, this reduction could result from two dif- 
ferent effects. It might be expected that doping will in- 
crease the interband scattering. Theory predicts that as 
the interband scattering is increased the large gap would 
decrease, while the small gap increases until they reach 
a common value and merge together. However, experi- 
mental studies (for example point contact spectroscopy)- 
have shown that the two gaps remain distinct, at least for 
low doping, and while doping appears to reduce the size 
of the large gap, the small gap remains constant. Alter- 
natively, the reduction of T c could follow from the effects 
of electron doping on the electronic structure. A recent 
theoretical study^ has concluded that the reduction of 
T c with increasing Al/C substitution can be explained 
by a reduction of the electron-phonon coupling constant 
which follows from a reduction in the density of states 
produced by electron doping. 

Although both Al and C doping produce a similar de- 
crease in T c , their effect on the upper critical held is 



very different. Al doping causes only a small change in 
H c2 and its anisotropy, however, C doping increases H C 2 
dramatically. 6 Similarly, C doping increases the residual 
resistivity at a much higher rate than Al doping. 6 - 7 The 
a band has most weight on the B plane, whereas the tt 
band has a large weight on both the B and Mg planes. 
It is therefore to be expected that replacing B with C 
should strongly increase the scattering rates on both the 
a and tt bands, but replacing Mg with Al mainly affects 
the 7T band. 

The de Haas-van Alphen (dHvA) effect is a powerful 
probe of the electronic structure of metals as it gives 
quantitative k-resolved information on the Fermi surface 
properties. In this paper, we present a dHvA effect study 
which quantifies the effect of the dopant on the electronic 
structure of Al doped MgB2 . 

Single crystal samples of Al and C doped MgB2 were 
prepared in Zurich using high pressure synthesis^ Pure 
MgB 2 crystals produced by this method have a T c of 
38.5 K, which is slightly (~0.5 K) lower than the best 
polycrystalline samples. 8 We measured two Al doped 
samples from batches AN215 and AN217 which both 
have T c ~ 33.6 K (see below). Clear dHvA signals 
were observed from both the Al doped samples for fields 
greater than 19 T. The Al content of the crystals was de- 
termined by measuring the c-axis lattice constants which 
have been shown 7 - to decrease linearly with increasing x 
(c = 3.513(2) - 0.28(2)x A). The two crystals were mea- 
sured to have c=3.4904(7)A and c=3.4920(10)A, giving 
x = 7.9 ±0.4% and x = 7.4 ±0.4% for AN215 and AN217 
respectively. Two C doped samples from batches AN314 
and AN284 which nominally have y = 3% and y = 4% 
with T c values of 35.7 K and 34.5 K respectively were also 
measured. No dHvA signals were observed from either 
sample for fields up to 33 T. This is consistent with a 
stronger increase in scattering rates in C doped samples, 
as indicated by their high residual resistances (see Fig.Q] 
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FIG. 1: Heat capacity vs temperature for Al doped MgB2 
samples AN215 and AN217. The data for AN217 have been 
shifted vertically for clarity. Inset: In-plane resistivity vs tem- 
perature, 1=C doped (AN314), 2=C doped (AN284), 3=A1 
doped (AN217), 4=pure (AN77). 



and RefsQ and0). 

Following our dHvA study we measured the heat ca- 
pacity of the two Al doped samples to probe for any 
homogeneous second superconducting phases. Heat ca- 
pacity is a bulk probe and hence is a useful complemen- 
tary probe to the dHvA measurements as dHvA signals 
will normally only be observed from the homogeneous 
bulk of the sample. The heat capacity c p was mea- 
sured as a function of temperature by an a.c. technique 9 
in fields of up to 7 T. The sample was placed on a 
flattened 12^m, chromel-constantan thermocouple and 
heated with a modulated light source. 

In Fig. □ we show Ac p = c p (B = 0) - c p (B" c =6 T) 
for both samples. As the samples are in the normal state 
for B^ c > 3.5 T, this procedure isolates that part of c p 
arising from the superconducting transition. The sam- 
ple labelled AN217 was the actual sample (mass=1.4/ig) 
used for the dHvA study, whereas sample labelled AN215 
(mass=70/ig) was part of the same crystal from which 
the much smaller dHvA sample was broken off. Sample 
AN217 has a sharp single transition, with T c =33.6 K (on- 
set), (10-90% width = 0.9 K), and dT c /dB^ c = -8.5±0.1 
K/T. Sample AN215 has more structure in its transition 
with a shoulder corresponding to approximately 25% of 
the total. T c =34.4 K, (onset) (10-90% width = 1.45 K) 
and dT c /dB^ c = -8.8 ± 0.1 K/T, although the onset of 
the main part of the transition is at approximately the 
same temperature as AN217. Measurements on a pure 
MgB 2 sample (batch AN189) gave T c =38.3 K, (onset) 
(10-90% width = 0.3 K) and dT c /dB^ c = -9.8 ± 0.1 
K/T. The smaller dT c /dB values in the Al doped sam- 
ples almost exactly scale with their reduced T c . 

Although our a.c. technique is very sensitive even for 
very small samples, it does not give accurate absolute 
values of c p and further for the smallest samples there 
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FIG. 2: Observed dHvA frequencies versus angle for the two 
Al doped samples AN215 (■,♦) and AN217 (0). Dashed lines 
are F(ff) for pure MgB2, solid lines are pure F(6) data scaled 
to fit the current data. 

is a large addenda contribution. For these reasons the 
units of c p are not quoted. For sample AN215 the jump 
in c p at T c is 15% of the total which is the same as that 
found for the pure MgB 2 samples. For AN217 the jump 
is ~ 3%, which results from the much larger addenda 
relative to the contribution from this tiny sample, rather 
than a reduced superconducting volume fraction. 

Quantum oscillations were measured using piezore- 
sistive cantilevers^ connected to a room temperature 
a.c. Wheatstone bridge. The torque is related to the 
change in bridge resistance by a factor sa 10 -10 Nmf2 -1 , 
with a noise level roughly corresponding to a torque of 
- 10- 14 Nm. 

The cantilever resistance was measured as the field 
was swept (typically at 0.8 T/min) between 20 T and 
33 T. Sweeps were made at fixed angle 8 as the sample 
was rotated from B\\c to (approximately) B\\a. The raw 
data were fast Fourier transformed to evaluate the fre- 
quencies F of any dHvA oscillations present. In Fig. [21 
we show the observed dHvA frequencies versus angle for 
the two Al doped samples. Comparing these results to 
those from pure MgB2, 11 ' 12 we find that the two higher 
amplitude signals originate from the minimal and maxi- 
mal extremal orbits of the smaller of the two quasi-two- 
dimensional a sheets of the Fermi surface^ with frequen- 
cies labelled previously F\ and F 2 respectively! 11 ! 12 The 
values of F\ (0 = 0) and F 2 (8 = 0) are the same in both Al 
doped crystals, and are ~ 10 — 20% smaller than the cor- 
responding values for pure MgB 2 (see Table HJ. Studies 
on many different samples of pure MgB 2 (from different 
sources) have remarkably reproducible dHvA frequencies 
(typically to within 1%) and so the smaller values found 
here are significant. The a sheets have an approximately 
circular in-plane cross-section and a cosine c-axis disper- 
sion and so we can calculate their volumes from the mea- 
sured extremal areas. We therefore estimate that there 
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are 16 ± 2% less holes in our Al doped samples compared 
to pure MgB 2 . 

Calculationsi2ii£*i& of the dHvA frequencies in pure 
MgB 2 are slightly different from the experimental val- 
ues. It was showiii2iiaii& that all the frequencies can be 
brought into agreement with theory by rigidly shifting 
both a bands down by ~80 meV and both tt bands up 
by ~60 meV relative to the Fermi level. For our Al doped 
samples we find that both a band frequencies can be ex- 
plained by a ~150 meV shift downward relative to the 
calculations for the pure material (or a ~70 meV shift rel- 
ative to the experimental results for the pure samples). 

To describe the effect of doping theoretically we have 
used the the virtual crystal approximation (VCA) to 
calculate the electronic band structure self-consistently 
on a very dense 47 x 47 x 35 k-mesh. The dHvA fre- 
quencies were then calculated as described previously^ 
Here the effect of doping with Al is simulated by re- 
placing the Mg atom with a virtual atom with charge 
Z = xZal — (1 — x)ZMg- In order to account for the 
difference between the dHvA areas for pure MgB2 and 
the calculations, we applied the same constant shift for 
pure MgB2 as described above to all doping concentra- 
tions. The VCA calculations give the following depen- 
dence of dHv A frequencies on doping x, F\ — 540— 1820x, 
F% = 1530 — 2050x. Hence, the observed frequencies 
(Fi and F%) correspond to a doping of x = 7.5 ± 1% 
and 8.4 ± 1% respectively. These values compare fa- 
vorably to the x values deduced from the c-axis lattice 
constant (x = 7.7 ± 0.4%). Close to x = the VCA 
calculations (including phonon renormalization)^ predict 
dT c /dx = —0.50 K/%, and hence for our doping we ex- 
pect a T c reduction of 4.0±0.3 K. The actual T c reduction 
was 3.9 K, and so we conclude that the VCA results ac- 
curately describe the effect of (relatively light) doping in 
MgB 2 . 

In sample AN215, for angles close to = 90°, a further 
frequency (-F4) is observed which originates from the hole- 
like 7r sheet. This frequency is also smaller than that of 
pure MgB2, however, the relatively large error precludes 
any estimate of the doping on this sheet. We note that 
no signals were observed from the electron- like it sheet, 
which in pure MgB2 produces the strongest dHvA signal 
F 3 (RefeE3l2). 

The quasiparticle effective mass m* on orbit F2 was 
determined by measuring the temperature dependence of 
the amplitude of the dHvA oscillations and fitting to the 
usual Lifshitz-Kosevich formula, il These measurements 
were made at = 10° ± 1° and 6 = 20° ± 1° for AN215 
and AN217 respectively. A reduction of 1.5% and 7.0% 
respectively was applied to give the mass at 9 = 0. Here 
we have used the fact that for this orbit m* scales accu- 
rately with the dHvA frequency. 18 Results for the two Al 
doped samples along with several undoped samples are 
shown in Table [n] The amplitude of the signals from 
Fi and F4 in these samples are too small to perform an 
accurate mass determination. 

As we do not expect there to be strong electron- 



TABLE I: Calculated (LDA)i& and measured dHvA frequen- 
cies [F° = F(6 = 0)] for pure±£ and Al doped MgB 2 . AE is 
the rigid band shift needed to bring the theoretical values in 
line with experiment. 





LDA 


Pure 




AN215 


AN217 


Orbit 


F° 


F° 


AE 


F° 


AE 


F° 




[T] 


[T] 


[meV] 


[T] 


[meV] 


[T] 


Fi 


730 


546±20 


85 


410±20 


148 




F, 


1756 


1533±20 


93 


1360±20 


147 


1360±20 


F :i 


2889 


2685±20 


-75 








F A 


458 


553±10 


-45 


480±40 


-10 





electron interactions in MgB2 we can use our value of 
m* along with the calculated bare band mass nib to es- 
timate the strength of the clcctron-phonon coupling con- 
stant A on this orbit, A = m*/ m & — !• Previously we have 
shown that for pure MgB2 the values of A calculated by 
this method are in good agreement with results of band 
structure calculations ^ILll A complication here is that 
nib varies slightly with the doping. In Table UTI we show 
the VCA calculated band masses for the pure and Al 
doped samples, along with the calculated A values. Al- 
though there is some variation between the values for the 
pure samples, A for the AN215 is 10 ± 5% smaller (the 
error for AN217 is too large for meaningful comparison). 

The theoretical electron-phonon mass enhancement on 
the a bands is the sum of intraband and interband parts, 
which for the undoped case are \ aa — 1.02 and X a7r = 
0.21 respectively, giving Af ot = 1.23 (note that AjT ot does 
not vary substantially within each a sheet and only varies 
by ~ 7% between the two a sheets) ML For £=7.7% the 
VCA calculation gives A^ ot = 1.16, which is ~ 5% smaller 
than for x — 0. This is in line with our observations. 

The field dependence of the dHvA amplitude r osc is 
proportional to B^RtRd- Evaluating the Dingle factor 
Rd = exp[— Trhkp / (eB£)], where kp is the orbitally aver- 
aged Fermi wavevector, allows us to estimate the mean- 
free-path £ on a particular orbit. In Fig.[3]we show Dingle 
plots [log(T osc / (B 2 Rt) versus l/B] for the frequency F 2 
in both samples at selected angles. The data for AN217 



TABLE II: Summary of m* and A values for F2 for several 
pure crystals (B and K from Ref ll2t) . AN77 (this work), and 
the two Al doped samples. Note that sample K contained iso- 
topically pure 10 B, whereas the others contain natural (mixed 
isotope) B. 



Sample 


m* 


nit 


A 


B 


0.65±0.01 


0.305 


1.12±0.03 


K 


0.61±0.01 


0.305 


1.00±0.03 


AN77 


0.61±0.02 


0.305 


0.99±0.06 


AN215 


0.57±0.01 


0.297 


0.92±0.03 


AN217 


0.59±0.03 


0.297 


0.97±0.1 
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FIG. 3: (a) Field dependence of the oscillatory torque (Fourier 
filtered to remove frequencies below lkT). AN215 data have 
been multiplied by 3 for clarity, (b) Raw oscillatory torque 
for fields close to the basal plane showing the contribution 
from F4 (c) Dingle plots for F2 for both samples. 

is exponential and is fitted by Rd = exp(— 165/S), giv- 
ing I — 270 A. This is approximately half of that found 
for pure samples of MgB2 (grown with natural mixed- 
isotope boron) The data for AN215 is markedly non- 
exponential and is consistent with there being a beat be- 
tween two different dHvA signals with frequencies differ- 
ing by ~40 T. As such behavior was not found for AN217 
or pure MgB 2 this indicates that there are two macro- 
scopic parts to AN215 which have different Al contents, 



in agreement with the specific heat measurements If 
two frequencies are included in the fit to the data then 
the scattering rate is found to be close to that of AN217. 

The small amplitude of F4 prevents us from determin- 
ing I accurately, however its relative amplitude suggests 
an approximate two-fold decrease in t. The signal from 
the hole- like it sheet (frequency F3) is normally large for 
fields oriented near the planefii*i£ and as we have not 
observed this orbit in either Al doped sample we con- 
clude that that the mean free path on this sheet must be 
reduced by a somewhat larger factor (> 3). 

In conclusion, we have measured dHvA oscillations 
in Al doped MgB 2 . Our results shows that the reduc- 
tion in size of the a sheets is in good agreement with 
band structure calculations using the VCA approxima- 
tion. The measured reduction in quasiparticle effective 
mass is also consistent with the reductions in band mass 
and electron phonon coupling predicted by the VCA cal- 
culations. These calculations also correctly explain the 
magnitude of the reduction of T c found in our samples. 
As expected Al doping increases the scattering more on 
the 7T bands than on the a bands. This excellent agree- 
ment between the dHvA results and theory can be seen 
as a direct confirmation of the effect of band filling on 
the superconducting properties which has been proposed 
theoretically. 4 

We thank I.I. Mazin for useful discussions and B.J. 
Pullum and G.M. Armstrong for technical assistance. 
This work was supported by the EPSRC (U.K.), and the 
NSF through grant number NSF-DMR-0084173. LB ac- 
knowledges support from the NHMFL in-house program. 



1 A. Y. Liu, I. I. Mazin, and J. Kortus, Phys. Rev. Lett. 87, 
087005 (2001). 

2 H. J. Choi, D. Roundy, H. Sun, M. L. Cohen, and S. G. 
Louie, Nature 418, 758 (2002). 

3 I. I. Mazin, O. K. Andersen, O. Jepsen, O. V. Dol- 
gov, J. Kortus, A. A. Golubov, A. B. Kuz'menko, and 
D. Van der marel, Phys. Rev. Lett. 89, 107002 (2002). 

4 J. Kortus, O. V. Dolgov, R. K. Kremer, and A. A. Golubov, 
Phys. Rev. Lett. 94, 027002 (2005). 

5 R. S. Gonnelli, D. Daghero, G. A. Ummarino, A. Calzolari, 
V. Dellarocca, V. A. Stepanov, S. M. Kazakov, J. Jun, and 
J. Karpinski, cond-mat/0407267 ((2004)). 

6 S. M. Kazakov, R. Puzniak, K. Rogacki, A. V. Mironov, 
N. D. Zhigadlo, J. Jun, C. Soltmann, B. Batlogg, and 
J. Karpinski, Phys. Rev. B 71, 024533 (2005). 

7 J. Karpinski, N. D. Zhigadlo, G. Schuck, S. M. Kazakov, 
B. Batlogg, K. Rogacki, R. Puzniak, J. Jun, E. Miiller, 
P. Wagli, R. Gonnelli, D. Daghero, G. A. Ummarino, and 
V. A. Stepanov, Phys. Rev. B 71, 174506 (2005). 

8 J. Karpinski, S. M. Kazakov, J. Jun, M. Angst, R. Puzniak, 
A. Wisniewski, and P. Bordet, Physica C 385, 42 (2003). 

9 A. Carrington, C. Marcenat, F. Bouquet, D. Colson, 
A. Bertinotti, J. F. Marucco, and J. Hammann, Phys. Rev. 
B 55, R8674 (1997). 

J. R. Cooper, A. Carrington, P. J. Meeson, E. A. Yelland, 



N. E. Hussey, L. Balicas, S. Tajima, S. Lee, S. M. Kazakov, 
and J. Karpinski, Physica C 385, 75 (2003). 

1 E. A. Yelland, J. R. Cooper, A. Carrington, N. E. Hussey, 
P. J. Meeson, S. Lee, A. Yamamoto, and S. Tajima, Phys. 
Rev. Lett. 88, 217002 (2002). 

2 A. Carrington, P. J. Meeson, J. R. Cooper, L. Bali- 
cas, N. E. Hussey, E. A. Yelland, S. Lee, A. Yamamoto, 
S. Tajima, S. M. Kazakov, and J. Karpinski, Phys. Rev. 
Lett. 91, 037003 (2003). 

3 J. Kortus, 1. 1. Mazin, K. D. Belashchenko, V. P. Antropov, 
and L. L. Boyer, Phys. Rev. Lett. 86, 4656 (2001). 

4 H. Rosner, J. M. An, W. E. Pickett, and S. L. Drechsler, 
Phys. Rev. B 66, 024521 (2002). 

5 H. Harima, Physica C 378, 18 (2002). 

6 I. I. Mazin and J. Kortus, Phys. Rev. B 65, 180510 (2002). 

7 D. Shoenberg, Magnetic Oscillations in Metals (Cam- 
bridge University Press, Cambridge, 1984). 

8 See EPAPS Document No. XXXX. 

9 Note that by scaling F2 with the decrease in the bulk T c , 
a 40 T difference in F2 corresponds to a difference in T c of 
~1 K in agreement with the c p data. 



